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Abstract 

Background: Bean anthracnose is caused by the fungus Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.- 
Scrib. Resistance to C lindemuthianum in common bean {Phaseolus vulgaris L.) generally follows a qualitative mode 
of inheritance. The pathogen shows extensive pathogenic variation and up to 20 anthracnose resistance loci 
(named Co-), conferring resistance to specific races, have been described. Anthracnose resistance has generally 
been investigated by analyzing a limited number of isolates or races in segregating populations. In this work, we 
analyzed the response against eleven C. lindemuthianum races in a recombinant inbred line (RIL) common bean 
population derived from the cross Xana x Cornell 49242 in which a saturated linkage map was previously 
developed. 

Results: A systematic genetic analysis was carried out to dissect the complex resistance segregations observed, 
which included contingency analyses, subpopulations and genetic mapping. Twenty two resistance genes were 
identified, some with a complementary mode of action. The Cornell 49242 genotype carries a complex cluster of 
resistance genes at the end of linkage group (LG) Pvl 1 corresponding to the previously described anthracnose 
resistance cluster Co-2. In this position, specific resistance genes to races 3, 6, 7, 19, 38, 39, 65, 357, 449 and 453 
were identified, with one of them showing a complementary mode of action. In addition, Cornell 49242 had an 
independent gene on LG Pv09 showing a complementary mode of action for resistance to race 453. Resistance 
genes in genotype Xana were located on three regions involving LGs PvOl, Pv02 and Pv04. All resistance genes 
identified in Xana showed a complementary mode of action, except for two controlling resistance to races 65 and 
73 located on LG PvOl, in the position of the previously described anthracnose resistance cluster Co-1. 

Conclusions: Results shown herein reveal a complex and specific interaction between bean and fungus genotypes 
leading to anthracnose resistance. Organization of specific resistance genes in clusters including resistance genes 
with different modes of action (dominant and complementary genes) was also confirmed. Finally, new locations for 

anthracnose resistance genes were identified in LG Pv09. 

Keywords: Common bean, Phaseolus vulgaris, Colletotrichum lindemuthianum, Anthracnose resistance inheritance. 
Genetic analyses. Genetic linkage map 
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Background 

Plants can recognize potential pathogens via two per- 
ception systems. One of them, named pathogen- or 
microbe-associated molecular patterns (PAMPs or MAMPs, 
respectively), detects conserved molecules associated 
with groups of pathogens through pattern recognition 
receptors leading to PAMP-triggered immunity. The other 
evolved to recognize pathogen virulence effectors, usually 
through intracellular resistance proteins (R proteins), 
causing effector-triggered immunity (ETI). ETI corre- 
sponds to what is classically referred to as gene-for- 
gene, vertical or race-specific resistance [1-3]. One of 
the first examples of race-specific resistance in plants 
was described by Barrus [4,5] in the interaction of Colleto- 
trichum lindemuthianum and common bean (Phaseolus 
vulgaris L.). 

Anthracnose, caused by the ascomycete fungus C. 
lindemuthianum (Sacc. & Magnus) Lams.- Scrib., can 
result in a devastating disease in common bean. Bean 
anthracnose has a worldwide distribution but it is par- 
ticularly problematic in temperate regions. The pathogen 
can attack all aerial parts of the bean plant and produces 
round black shrunken lesions containing flesh colored 
spores on leaves, stem, pods and seeds. The attack of 
this fungus can result in premature defoliation, prema- 
ture fall of flowers and pods, seed deterioration and, in 
extreme cases, plant death. Infected seeds are the major 
means of dispersal of the pathogen [6]. The pathogen 
shows extensive pathogenic variation, with at least 100 
pathogenic variants or races reported among isolates of 
C. lindemuthianum collected worldwide [7-11] using a 
set of 12 differential cultivars and a standardized method 
to name the races [12]. 

Resistance to C. lindemuthianum generally follows a 
qualitative mode of inheritance where resistant and sus- 
ceptible reactions are clearly differentiated. Identification 
of anthracnose resistance genes by classical genetics is 
based on the interpretation of results obtained from F2 
segregating populations derived from two types of cros- 
ses: R X S or R X R (R is resistant and S is susceptible). 
Results observed in R x S crosses are used to infer the 
number and mode of action of genes conferring resist- 
ance to C. lindemuthianum, while those for R x R cros- 
ses are used to identify the specific genes involved in 
the reaction against this pathogen (allelism tests). Since 
the first anthracnose resistance gene was reported [13], 
many genetic analyses have been conducted to study an- 
thracnose resistance inheritance in different bean geno- 
types. Up to 20 anthracnose resistance loci conferring 
resistance to specific races (designated as Co- followed 
by a number or a letter) have been described in common 
bean: Co-1 to Co- 7, co-8, Co-9 to Co-14 and Co-u to 
Co-z [14]. All anthracnose resistance genes described, 
except for co-8, exhibit complete dominance where the 



dominant allele conditions the resistance reaction. A 
complementary mode of action between two independ- 
ent resistance genes has also been described using F2 or 
F2:3 segregating populations, being necessary the pres- 
ence of both dominant alleles for expression of resist- 
ance [15-17]. Many genetic analyses assumed that the 
resistance to different races in a bean genotype is con- 
ferred by the same gene. Based on this hypothesis, most 
classical studies considered that different resistance 
spectra in genotypes were due to different alleles of the 
same gene. So, different alleles were described for genes 
Co-1, Co-3 and Co-4 [14]. 

Most identified anthracnose resistance genes were lo- 
cated in the genetic map of common bean: genes Co-1, 
Co-x and Co-w were mapped on linkage group (LG) 
PvOl [18,19]; Co-u was located on Pv02 [18]; Co-13 on 
Pv03 [20]; Co-2 on Pvll [21]; Co-3, Co-9, Co-y, Co-z and 
Co-10 on Pv04 [19,22-26]; Co-4 on Pv08 [19,27]; and 
Co-5, Co-6 and Co-v on Pv07 [28,29]. Mapping of genes 
conferring resistance to several specific races revealed 
that some of these Co- genes were organized in clusters 
of race-specific resistance genes. Tight linkage associa- 
tions of many Co- resistance genes have been well es- 
tablished on LGs Pv04, Pv07 and Pvll [15,19,26,28] 
corresponding to the named clusters Co-3, Co-5 and 
Co-2, respectively. At a molecular level, the majority of 
plant R genes cloned so far encode proteins with two 
specific domains: nucleotide-binding-sites (NBS) and a 
leucine-rich repeat (LRR) domain. Genes encoding R 
proteins were found in tandem on chromosome regions 
corresponding to the Co-3 and Co-2 clusters observed 
in genetic analyses [30-33]. 

The anthracnose resistance system in common bean 
has been classically investigated by analyzing a limited 
number of isolates or races in different segregating pop- 
ulations. In the present study, the response against 11 C. 
lindemuthianum races was analyzed in a recombinant 
inbred line (RIL) population derived from the cross 
Xana x Cornell 49242 in which a saturated linkage map 
was previously developed [34,35]. The present study 
aims to add understanding concerning the organization 
and interaction between anthracnose resistance genes, 
and reveals a complex P. vulgaris-C. lindemuthianum 
interaction. 

Results 

Segregations for the eleven C. lindemuthianum races 

Parental line Xana was susceptible to races 6, 38, 39, 357 
and 453, and resistant to the remaining six races (3, 7, 
19, 65, 73 and 449). Parental line Cornell 49242 was sus- 
ceptible to race 73 and resistant to the remaining ten 
races. Table 1 shows the observed segregations for re- 
sistance to each race in the XC RIL population (see 
Additional file 1 for more detail of the segregation ratios 
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Table 1 Observed segregations for resistance to eleven 
Colletotrichum lindemuthianum races in the XC RIL 
population 

Parental phenotype XC RIL population 



Cornell Observed Segregation ratio 



Race 


Xana 


49242 


R 


S 


Tested (R: S) 




P 


6 


S 


R 


47 


38 


1: 1 


0.95 


0.33 


38 


S 


R 


48 


42 


1: 1 


0.40 


0.53 


39 


s 


R 


40 


39 


1: 1 


0.01 


0.91 


357 


s 


R 


37 


42 


1: 1 


0.32 


0.57 


73 


R 


S 


37 


40 


1: 1 


0.12 


0.73 


65 


R 


R 


82 


20 


3: 1 


1.58 


021 


3 


R 


R 


67 


35 


5: 3 


0.44 


0.51 


7 


R 


R 


65 


35 


5: 3 


0.27 


0.61 


19 


R 


R 


71 


35 


5: 3 


0.91 


0.34 


449 


R 


R 


64 


39 


5: 3 


0.01 


0.94 


453 


S 


R 


50 


47 


1: 1 


0.09 


0.76 



The segregation ratio considered and tile chi-square goodness-of-fit test is 
indicated in each case. R, resistant; S, susceptible. 



expected). Segregations for resistance to races 6, 38, 39, 
357, 73 and 453 conformed with the 1:1 resistant:sus- 
ceptible (R:S) ratio, expected for one resistance gene or 
for three independent resistance genes, complementary 
two-by-two. Segregation for resistance to race 65 con- 
formed with the 3:1 R:S ratio, expected for two inde- 
pendent genes. Segregations for resistance to races 3, 7, 
19 and 449 conformed with the 5:3 R:S ratio, expec- 
ted when three independent genes, two with a comple- 
mentary mode of action, are involved in the resistance 
response. 



Genetic analyses of resistance to races 6, 38, 39 and 357 

Segregation for resistance to races 6, 38, 39 and 357 
[Xana (S) x Cornell 49242 (R)] conformed with the 1:1 
R:S ratio (p = 0.33, p = 0.53, p = 0.91, p = 0.57; Table 1). 
Contingency chi-square tests corresponding to the joint 
segregations for each of the resistances with markers 
tagging the six main anthracnose resistance clusters de- 
viated significantly from what was expected by random 
in the case of markers SQ4 and SCAreoli, that tagged 
the Co-2 resistance cluster on LG Pvll (Table 2). 

To determine if a single gene or different genes of 
the same cluster were involved in the response, co- 
segregations for resistance to races 6, 38, 39 and 357 
were considered. Co-segregation for resistance to the 
four races was observed in 69 RILs (35 RILs were re- 
sistant to the four races, R*^ R^** R^'' R^^'; and 34 were 
susceptible, S^ S^* S^^ S^^^), and in three cases evidence of 
recombination was found (one RIL with each of the fol- 
lowing haplotypes: R*^ S^« R^*^ S^^^, R'^ R^« R^'^ S^*^^ and R*- 
g38 g39 g357-j Accordingly, it can be concluded that four 
different closely linked genes at the Co-2 resistance cluster 
in Cornell 49242 determined specific resistance to races 6, 
38, 39 and 357 (genes Co-^'^, and Co- 

2^57 -Cy gjj^^g ^jj evidence suggests a monogenic control of 
resistance to each of the four races, these resistance genes 
were directly included in the genetic map (Figure 1). 

Genetic analyses of resistance to race 73 

Resistance to race 73 [cross Xana (R) x Cornell 49242 
(S)] conformed with the 1:1 R:S segregation ratio (p = 0.73; 
Table 1). The contingency chi-square values (Table 2) de- 
viated significantly from the expectation of random seg- 
regation when compared with markers CV542014 and 
OFIO530, tagging the Co-1 resistance cluster. This result 



Table 2 Results of contingency chi-square tests conducted between the joint segregations for each one of the resistance 
to eleven races and eleven loci tagging the main anthracnose resistance clusters identified in common bean 









Race6 


Race38 


Race39 


Race357 


Race73 


Race65 


Race3 


Race7 


Race19 


Race449 


Race453 


Marker 


LG 


Cluster 


X^-P 


X^-P 


X^-P 


X^-P 


X^-P 


x'-P 


X^-P 


X^-P 


x'-P 


X^-P 


x'-P 


CV542014 


PvOl 


Co-1 


0.10 ns 


1.76 ns 


1.96 ns 


3.15 ns 


44.62* 


12.87* 


1.61 ns 


2.30 ns 


2.20 ns 


3.16 ns 


2.03 ns 


OFIO530 


PvOI 


Co-1 


0.07 ns 


0.54 ns 


038 ns 


0.88 ns 


47.69* 


14.95* 


0.48 ns 


1 40 ns 


0.60 ns 


1.36 ns 


1.25 ns 


/-gene 


Pv02 


Co-u 


1 .22 ns 


0.62 ns 


1.58 ns 


0.90 ns 


0.04 ns 


0.01 ns 


0.49 ns 


0.23 ns 


0.04 ns 


0.07 ns 


0.51 ns 


SWT 3 


Pv02 


Co-u 


3.57 ns 


1.68 ns 


2.84 ns 


0.66 ns 


0.11 ns 


0.12 ns 


0.04 ns 


0.19 ns 


0.42 ns 


0.07 ns 


0.00 ns 


254-G15F550 


Pv04 


Co-3 


0.32 ns 


0.00 ns 


0.06 ns 


0.06 ns 


0.65 ns 


0.20 ns 


15.62* 


13.11* 


1411 


11.85 ns 


1 7.02* 


SWT 2 


Pv04 


Co-3 


041 ns 


0.05 ns 


0.31 ns 


0.63 ns 


021 ns 


1 .83 ns 


1 7.99* 


19.11* 


16.58 


13.91* 


26.24* 


Pha 


Pv07 


Co-5 


0.10 ns 


0.00 ns 


0.59 ns 


0.17 ns 


1.19 ns 


0.72 ns 


0.17 ns 


0.00 ns 


0.09 ns 


0.10 ns 


0.31 ns 


SZ4b 


Pv07 


Co-5 


0.01 ns 


0.17 ns 


0.06 ns 


0.10 ns 


2.06 ns 


1.11 ns 


0.03 ns 


0.05 ns 


0.08 ns 


0.00 ns 


3.14 ns 


SBB14 


Pv08 


Co-4 


0.00 ns 


0.04 ns 


0.05 ns 


0.13 ns 


1.47 ns 


1.11 ns 


0.12 ns 


0.07 ns 


0.35 ns 


0.02 ns 


0.03 ns 


SQ4 


Pvll 


Co-2 


52.60* 


57.96* 


45.06* 


51.19* 


2.11 ns 


1 1 43 ns 


32.54* 


32.75* 


38.09* 


35.56* 


12.45* 


SCAreoli 


Pvll 


Co-2 


56.44* 


62.82* 


51.11* 


57.74* 


2.31 ns 


13.33 


35.24* 


39.13* 


41.90* 


38.89* 


14.54* 



LG = linkage group; ns = not significant, * = significant considering a Bonferroni-correted significance threshold level of 0.0005 for multiple comparisons. 
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Figure 1 Bean linkage groups in which anthracnose resistance genes were directly or indirectly located. Five race-specific resistance 
genes were mapped on LGs PvOI, in a position corresponding to that of Co-1 resistance cluster (gene CoV^'^), and Pvl 1, in a position corresponding 
to that of Co-2 cluster (genes and Ten resistance genes were located on LGs Pv02 (CoPvOlc''''' , CoPv02c'''', 

CoPvOlc""" CoPvOci"^''), Pv04, in a position corresponding to that of Co-3 cluster [Co-3(^''', Co-Bc^'", Co-Sc'''''', Co-Sc"^'", and Co3c''^^~''), and Pv09 
{CoPv09c''^^''~). Resistance genes are named by using its location in the genetic map (LG or Co-anthracnose resistance cluster), name of the isolate or 
race (in superscript), followed by the bean genotype in which the resistance gene was identified (X, Xana; C, Cornell 49242). Genes with a complementary 
mode of action are indicated using the letter 'c' after the name of the gene.. Map distances, on the left, are expressed in centlMorgans, estimated using 
the Kosambi mapping function. M_ , AFLP marker loci; Sp, seed protein marker loci; conting, fZ_^ BM_ and PV, microsatellite marker loci; S_ and 2S4_, SCAR 
marker loci; Fin, locus controlling indeterminate versus determinate growth habit; /, gene controlling resistance against bean common mosaic virus; 
0_. RAPD markers; IND_, InDel markers. 



indicates that resistance to race 73 in Xana may be medi- 
ated by one resistance gene located at the Co-1 
cluster. Genetic mapping confirmed the position of this 
gene at the end of LG PvOI, in a position corresponding 
to cluster Co-1 (Figure 1). 

Genetic analyses of resistance to race 65 

Segregation of resistance to race 65 [cross Xana (R) x 
Cornell 49242 (R)] suggests the presence of two inde- 
pendent dominant genes (Table 1). The chi-square values 
deviated significantiy from the expectation of random seg- 
regation when resistance to race 65 was compared with 
molecular markers that tagged the Co-1 and Co-2 regions 
(Table 2). This result suggests the localization of one re- 
sistance gene at the Co-1 cluster, and a second gene 
at the Co-2 cluster. To confirm this model, segrega- 
tion for resistance to race 65 was analyzed in six subpopu- 
lations established from the total XC RIL population (see 



Additional file 2): X-Co-1 and C-Co-1, formed by 30 and 
47 RILs showing, respectively, the Xana and Cornell49242 
alleles for the markers CV542014 and OFIO530. X-Co-2 
and C-Co-2, formed by 42 and 41 RILs showing, respect- 
ively, the Xana and Cornell49242 alleles for the markers 
SQ4 and SCAReoli. X-Co-3 and C-Co-3, formed by 37 
and 41 lines, respectively, the Xana and Cornell49242 al- 
leles for markers 254-G15F550 and SW12. In the X-Co2 
subpopulation, the resistance to race 65 segregated ac- 
cording to a 1:1 R:S ratio, expected in the case of one re- 
sistance gene. A linkage analysis carried out within this 
subpopulation revealed a close linkage between this segre- 
gating gene and Co-f^'^ (RF = 0.00; LOD = 10.54). The 
localization of a resistance gene for race 65 derived from 
Xana at the Co-1 cluster was confirmed in the X-Col sub- 
population, in which the Co-1 region from Xana was 
fixed; the 27 RILs of this subpopulation were resistant to 
race 65. The opposite subpopulation, C-Col, showed a 
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segregation for resistance to race 65 according to the ratio 
expected for one dominant gene (i.e. 1:1 R:S). This resist- 
ance gene exhibited a close linkage with (RF = 
0.04; LOD = 7.68). The localization of a resistance gene 
against race 65 derived from Cornell 49242 at the Co-2 
cluster was confirmed in the C-Co2 subpopulation, in 
which the majority of lines showed a resistant phenotype 
against race 65 (all except for two lines, probably due to 
recombination). Independence between resistance against 
race 65 and the Co-3 chromosome region was deduced 
from subpopulations X-Co3 and C-Co3 - in both these 
subpopulations there was a 3:1 R:S ratio for resistance to 
race 65, as for the total XC population. 

It can be concluded that resistance to race 65 was con- 
ferred by one gene {Co-f^'^) in Xana located at the Co-1 
cluster on PvOl and by a gene in Cornell 49242 

located at the Co-2 cluster on Pvll (see Figure 1). 

Genetic analyses of resistance to races 3, 7, 19 and 449 

Segregation of resistance to races 3, 7, 19 and 449 [Xana 
(R) X Cornell (R)] conformed with the 5:3 R:S ratio {p = 
0.51, p = 0.61, p = 0.34, p = 0.94 respectively; Table 1) ex- 
pected for three independent genes, two of them having 
a complementary mode of action. Resistance to four 
races showed close co-segregation. A total of 96 lines 
showed the parental genotypes (63, R^ R'^ R^^ R**'' and 
33, S^S'S^^S**^), while only four lines revealing evidences 
of recombination (2, R^' R^" R*^"; and 2, R^ R' R^" S**'). 

Concerning resistance to race 3, the contingency chi- 
square values deviated significantly from the expectation 
of random segregation when resistance to this race was 
compared with molecular markers tagging the Co-3 and 
Co-2 chromosome regions (Table 2). This result can be 
interpreted as one of the resistance genes being located 
at the Co-3 cluster and a second gene located at the Co- 
2 cluster. Concerning the third complementary gene, it 
segregated independently from the molecular markers 
tagging the Co-1, Co-5, Co-4 and Co-u regions (Table 2), 
so contingency chi-square tests were conducted between 
resistance to race 3 and the remaining 368 loci included 
in the XC genetic map. There was significant deviation, 
with a block of three markers closely mapped on LG Pv02 
(Figure 1): IND_2_403966 (Cont. x^ = 4.24, p = 0.04); 
IND_2_404411 (Cont. x^ = 4.51, p = 0.03) and MctaEta^*^*^ 
(Cont. = 5.24, p = 0.02). 

Analyses within subpopulations (see Additional file 2) 
were consistent with this scenario. Within subpopula- 
tions involving the Co-1 cluster (X-Col and C-Col) re- 
sistance to race 3 fitted a 5:3 R:S ratio, as for the total 
XC RIL population, suggesting that the Co-1 region was 
not involved in genetic control against race 3. However, 
a change in the segregation ratio was observed within 
the two subpopulations involving the Co-3 cluster. Good 
fits to 3:1 and 1:1 R:S ratios were observed in the X-Co-3 



and C-Co-3 subpopulations, respectively. This is the ex- 
pected situation if one complementary gene proceeding 
from Xana is located at the Co-3 cluster. Concerning the 
Co-2 cluster, most lines included in the C-Co-2 subpopu- 
lation were resistant to race 3 (except one susceptible line, 
probably due to recombination); and the opposite X-Co2 
subpopulation fitted a 1:3 R:S ratio, expected for two inde- 
pendent and complementary genes (see Additional fUe 1 
for more detail). Both results are expected if the resistance 
gene derived from Cornell 49242 without a complemen- 
tary mode of action was located at the Co-2 cluster. Two 
subpopulations involving LG Pv02 were also considered 
for race 3 using markers IND_2_403966 and MctaEta^*^*^: 
the X-Pv02 subpopulation including 26 RILs showing the 
Xana allele for both markers, and the C-Pv02 subpopula- 
tion including 24 RILs showing the Cornell 49242 alleles. 
In the X-Pv02 subpopulation resistance to race 3 fitted a 
3:1 R:S ratio (21:5 R:S; x^ = 0.46, p = 0.50), while the C- 
Pv02 subpopulation fitted a 1:1 R:S ratio (12:12 R:S; x^ = 
0.00, p = 1.00). This situation supports the localization of a 
complementary resistance gene from Xana against race 3 
in this position of LG Pv02. 

In summary, results indicate involvement of three resist- 
ance genes against race 3 (Figure 1): one gene (€0-2^'*") 
from Cornell 49242 located at the Co-2 cluster, and two 
independent resistance genes from Xana with a comple- 
mentary mode of action, one (Co-3c^'^) located at the Co- 
3 cluster on Pv04 and a second {CoPv02c^'^) on LG Pv02. 
The same model deduced for resistance to race 3 is valid 
for resistance to races 7, 19 and 449, all genes showing 
close co-segregation. 

Genetic analyses of resistance to race 453 

Resistance to race 453 [cross Xana (S) x Cornell 49242 
(R)] fitted a 1:1 R:S segregation ratio {p = 0.76; Table 1). 
In this case, the chi-square values (Table 2) deviated sig- 
nificantly from the expectation of random segregation 
when compared with molecular markers that tagged two 
different chromosome regions: SW12 and OFlOnoo (Co-3 
region) and SQ4 and SCAreoli (Co-2 region). This finding 
provides evidence concerning the presence of more than 
one gene in the resistance response against race 453. Gi- 
ven the observed ratio, the most probable explanation is 
that resistance to race 453 was controlled by three inde- 
pendent resistance genes having a complementary mode 
of action two-by-two (see Additional file 1 for more de- 
tail). According to the results (Table 2), one of these genes 
would be located at the Co-2 region and another probably 
located at the Co-3 region. The third complementary gene 
segregated independently from molecular markers tagging 
the Co-1, Co-5, Co-4 and Co-u regions, so contingency 
chi-square tests between the resistance to race 453 and 
the remaining 368 loci included in the XC genetic map 
were conducted. A significant deviation was observed with 



Campa et al. BMC Plant Biology 2014, 14:115 
http://www.biomedcentral.com/1471-2229/14/115 



Page 6 of 1 2 



a block of 11 markers mapped on LG Pv09, from marker 
BM202 to McatEag^^* (see Figure 1 and Table 3). 

To confirm this model, segregation for resistance to 
race 453 was analyzed within the subpopulations estab- 
lished (Additional file 2). Segregation for resistance in 
the X-Col and C-Col subpopulations fitted a 1:1 R:S ra- 
tio, as for the total RIL population, suggesting that this 
region was not involved in genetic control of the reac- 
tion to race 453. In contrast, a change in the segregation 
ratio was observed in the remaining subpopulations. In 
the X-Co2 and C-Co-3 subpopulations resistance to race 
453 fitted a 1:3 R:S ratio, expected in the case of two 
complementary and independent genes. This result sug- 
gests that both subpopulations lacked one of the three 
complementary resistance genes. The X-Co2 subpopula- 
tion lacked a complementary resistance gene from Cor- 
nell 49242 located on the Co-2 region while the C-Co3 
subpopulation lacked a complementary resistance gene 
from Xana, located at the Co-3 region. The observed 
segregation within the opposite subpopulations, C-Co2 
and X-Co3, conformed with the 3:1 R:S ratio, thus being 
consistent with this scenario. If one of the three comple- 
mentary genes was fixed in a subpopulation, then a 3:1 
R:S ratio would be expected. 



Table 3 Genotype in the parental Xana and the 
recombinant inbred lines XC30, and XC88 for a block of 
markers in LGs Pv04, Pv09 and Pvl 1 of lines selected in 
order to verify the location and mode of action of the 
resistance genes deduced against race 453 

Region 
or cluster 



Race 453 Genotype^ 



Co3 



Co2 



LG 


Marker loci 


cM 


Cont.x^ 


P 


Xana 


XC30 


XC88 


Pv04 


SW12 




26.24 




A 


B 


B 


Pv04 


FZ-E9b 


2.7 


18.36 




A 


A'' 


B 


Pv04 


254-G15F^^° 


64 


8.52 




A 


B 


A'^ 


Pv04 


254-G15F''™ 


5.1 


16.62 




A 


B 


B 


Pv04 


254-G15F"" 


141 


17.02 




A 


B 


B 


Pv09 


IND_9_280580 




4.50 




A 


A 


B 


Pv09 


IND_9_291822 


1.1 


8.51 




A 


A 


B 


Pv09 


IND_9_291976 


0 


8.17 




A 


A 


B 


Pvll 






14.43 




A 


B 


A 


Pvll 


€0-2'"-^ 


0.6 


15.62 




A 


B 


A 


Pvll 




0.6 


1827 




A 


B 


A 


Pvll 




14 


19.49 




A 


B 


A 


Pvll 


MctgEat^**' 


2.3 


16.44 




A 


B 


A 


Pvll 


SQ4 


5.1 


12.45 




A 


B 


A 


Pvll 


SCAreoli 


0.7 


14.54 




A 


B 


A 



LG = linkage group. cM, distance between loci in centimorgans. 
* = significant. 

= genotype of the parental line Xana; B = genotype of the parental 
line Cornell49242. 
'^recombinant genotypes. 



For the third resistance gene against race 453, with 
location estimated at LG Pv09 based on contingency 
chi-square tests (Table 3), two additional subpopulations 
were considered using markers IND_9_280580 and 
ATA217. In the X-Pv09 subpopulation, formed by 46 
RILs showing the Xana allele for both markers, resist- 
ance to race 453 fitted a 3:1 R:S ratio (16:9 R:S, x^s^ = 
1.61, p = 0.20). In the C-Pv09 subpopulation, formed by 
30 RILs showing the Cornell 49242 allele, resistance to 
race 453 fitted a 1:3 R:S ratio (15:24 R:S, x^ls = 3.77, p = 
0.05). This result is in agreement with the localization 
of a third complementary resistance gene from Cornell 
49242 at LG Pv09. 

In summary, results indicate the involvement of three 
complementary genes, complementary two-by-two in the 
genetic control of response to race 453 (Figure 1): one 
gene (CoSc*^^"'') from Xana located on the Co-3 region, 
one (Co2c*'^^''^) from Cornell 49242 located on the Co-2 
region and a third gene (CoPv09c'^^^''^) from Cornell 
49242 located on LG Pv09. 

Genetic dissection of resistance to race 453 

Genetic dissection was performed to verify the complex 
system of resistance observed against race 453. Three 
susceptible lines were selected according to their paren- 
tal genotype (Xana or Cornell 49242) for markers lo- 
cated in LGs Pv04, Pv09 and Pvll (Table 3), in which 
were tentatively localized the three complementary re- 
sistance genes. Parental line Xana was susceptible to 
race 453 and according to this model it carried the com- 
plementary gene Co-3c^^^'^ and lacked Co-2c*^^''" and 
CoPv09c^^^ '^ . RIL XC30 showed the Cornell 49242 
genotype for markers of the Co-3 cluster, and the Xana 
genotype for markers that tagged the other two regions, 
the Co-2 cluster and the Pv09 region. Based on this, 
susceptible line XC30 should have the Co-2c*^^ gene 
and lack Co-3c*'^^'^ and CoPv09c^^^''^ . Finally, RIL XC88 
showed the Cornell genotype for markers tagging the 
Pv09 region and the Xana genotype for markers of clus- 
ters Co-2 and Co-3. Thus, it should have CoPv09c*^^'^ 
and lack Co-3c*^^'^ and Co-2c'''^^''". 

A total of 13 Fi seedlings were obtained from the cross 
between the susceptible lines Xana (Co-3c ' ) x RIL 
XC30 {Co-2c*^^''~') and all were resistant to race 453, in- 
dicating a complementary mode of inheritance. A total 
of 16 Fi seedlings were obtained from the cross Xana 
(Co-3c*^^'^) X RIL XC88 (CoPv09c*^^^'^) and all were re- 
sistant to race 453. Finally, from the cross RIL XC88 
{CoPv09c^^^''^) X RIL XC30 (Co-2c*^^''=) a total of 23 seed- 
lings were obtained, all resistant to race 453. In all cases, 
cross authenticity was verified using codominant molecular 
markers. This result confirmed the model of three genes, 
complementary two-by-two, involved in the genetic control 
of resistance to race 453 in the XC RIL population. 
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Discussion 

In the present study, the inheritance of resistance to 11 
C. lindemuthianum isolates classified as different races 
was analyzed in a RIL population derived from the cross 
Xana x Cornell 49242. In classical genetic analysis, the 
identification of a resistance gene is based its relation- 
ship of independence or linkage with other genes previ- 
ously described, through allelism tests. Molecular markers 
can help in determining the identity of a resistance gene. 
Molecular markers linked to the main anthracnose resis- 
tance genes have been reported [14] and can be used to 
establish the identity of a resistance gene by means of link- 
age analyses. The use of saturated genetic linkage maps al- 
lows more accurate location of a resistance gene by direct 
mapping, and also the identification of new resistance loci, 
independent of those previously described. Nevertheless, 
the use of linkage maps for direct mapping is limited by 
the number of genes involved in the resistance, and/or by 
possible epistatic interaction between them. This scenario 
occurred for the anthracnose resistance genetic analyses 
conducted on the XC RIL population. Only reactions to 
races 6, 38, 39, 357 and 73 exhibited a monogenic se- 
gregation and the respective resistance genes involved 
were directly mapped. Resistance to the remaining six 
races showed complex segregations, involving several 
Mendelian genes with different modes of action. In 
this case, genetic analyses supported by contingency chi- 
square tests and subpopulation analyses were proposed as 
an alternative solution for the analyses of complex segre- 
gations involving more than one gene. Genome-wide ana- 
lysis using contingency chi-square tests were conducted to 
identify the regions of the genetic map showing a signifi- 
cant association with the reaction to a specific race. Sig- 
nificant associations between the response to a specific 
race and a mapped locus suggest a genetic linkage rela- 
tionship between them, particularly if significant devia- 
tions are caused by an excess of the parental classes. 
Subpopulation analysis allows simplification of segregation 
and testing the involvement of a specific region in the 
control of the resistance. However, subpopulation size is 
reduced in respect to that of the original population, so 
larger original populations are required. Finally, genetic 
dissection was performed in specific cases to check the 
position of a resistance gene and its mode of action. Ne- 
vertheless, genetic dissection requires the development of 
new segregating populations from genotypes showing a 
specific combination of alleles, and these genotypes are 
not always available within the original segregating po- 
pulation due to recombination events. This systematic 
genetic analysis has allowed the drawing of conclusions 
concerning genetic control of anthracnose resistance in 
parental lines Xana and Cornell 49242. 

Results indicate that Cornell 49242 carried a complex 
cluster of race-specific resistance genes at the end of LG 



Pvll, corresponding to the Co-2 cluster. Four resistance 
genes and were 

directly mapped on this cluster, closely linked to markers 
SQ4 and SCAreoli. Recombinant lines in the response to 
the four races were detected, indicating that different re- 
sistance genes were present. From the analyses of sub- 
populations, the presence of another six resistance genes 
at this position was indirectly deduced: Co-2^''~', Co-2'''~', 

ol9-C ^ r,65-C J r' n449-C 11 ^ o 453-C 

Co-2 , Co-2 and Co-2 , as well as Co-2c 
with a complementary mode of action. The dominant 
resistance gene Are (re-named Co-2) was reported in 
Cornell protecting against races alpha, beta and 

gamma [36]. Resistance genes showing a complementary 
mode of action have been also suggested in this geno- 
type from the observed segregations in F2 populations 
[17]. Gene Are was introgressed into cultivar Ms8E02 
from Cornell 49242 [21,37]. Using a BCiFj population 
obtained from the backcross (Ms8E02/*Corel), gene Are 
was mapped at the end of LG PI (corresponding to LG 
Pvll) closely linked to the SCAR marker SCAreoli [37]. 
Results of the present study indicated that the original 
gene Are described in Cornell 49242 was in fact a cluster 
of linked resistance genes with each one conditioning re- 
sistance to different races of C. lindemuthianum. A clus- 
ter organization including closely linked genes against 
races 6, 31, 38, 39 and 65, was also identified in the 
A252 bean genotype in the same genetic position using 
the F2:3 population derived from the cross Andecha x 
A252 [19]. 

Results also indicated that genotype Cornell 49242 car- 
ried two complementary resistance genes to race 453, 
one [Co-2c^'^^^'^) located in the described cluster Co-2, 
and another {(ZoPv09c'^^'^) in LG Pv09 between mar- 
kers BM202 and ATA217 - this is the first anthracnose 
resistance gene located in this LG. Interestingly, a quan- 
titative trail locus (QTL) involved in resistance against 
Pythium ultimum [38] was localized to this same relative 
position. Evidence supplied by in silico analysis in the se- 
quenced bean genotype G 19833 (www.phytozome.net) 
was consistent with the location of anthracnose resist- 
ance genes in this position. At least four genes encoding 
NBS-LRR proteins and four genes encoding protein ki- 
nases, all typical R proteins, were annotated between 
markers BM202 and ATA217 (see Additional file 3). 

Genetic control of anthracnose resistance in the culti- 
var Xana had not been previously analyzed. Xana showed 
a resistant reaction against six of the eleven isolates ana- 
lyzed. Genetic analysis suggested that regions on LGs 
PvOl, Pv02, Pv04 and Pvll were implicated in the re- 
sponse. Resistance to race 73 showed monogenic inherit- 
ance and the corresponding locus {Co-f^'^) was directly 
mapped on the end of LG PvOl. Subpopulation analysis 
allowed the deduction of the presence of one additional 
resistance gene {Co-f'^"^) in this relative position. Co-f'^'^ 
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probably corresponds to the resistance gene against race 
65 mapped in this same relative position from the geno- 
type Andecha [19], as Andecha is one of the parents from 
which Xana was obtained. The Co-i^^'^ gene is closely 
linked to markers OFIO530, TGAl.l and CV542014. The 
OFIO530 fragment was previously linked to the anthrac- 
nose resistance locus Co-1 in the ¥2-3 population derived 
from the cross between the near-isogenic lines N85006 S 
and N85006 R [39] and later mapped to PvOl [18,19,40]. 
A gene conferring resistance to race 73 in AND277 was 
also located on PvOl using a F2 population derived from 
the cross AND277 x Ruda [41]. This resistance gene in 
AND277, that corresponds to Co-1, is closely linked to 
markers CV542014 and TGAl.l. In the genotype Kaboon, 
Co-1 was also shown to have a cluster organization, in- 
cluding three resistance genes against races 31, 81 and 
1545 [15]. In conclusion, different genetic evidence sup- 
ported the presence of a cluster, including race-specific re- 
sistance genes to anthracnose in the relative position in 
which the Co-1 gene was mapped (named cluster Co-1). 

Resistance genes derived from Xana against races 3, 7, 
19 and 449 were indirectly located on LGs Pv02 and 
Pv04, all with a complementary mode of action. A 
gene conferring specific resistance to C. lindemuthianum 
strains E4 and E42b (showing complete co-segregation) 
was mapped to the end of LG Pv02 in the Mesoamerican 
genotype BAT93 using the RIL population BAT93 x 
JALOEEP558. This gene, named Co-u, was located in the 
vicinity of the / locus [18], a complex resistance clus- 
ter effective against potyviruses [42]. The correspond- 
ence between Co-u and the resistance cluster identified in 
Xana is not clear, because the / locus and the SW13 
marker (linked to gene I) are included in the XC genetic 
map, and segregated independently from resistance genes 
CoPv02c^^, CoPv02c'^, CoPv02c'^'^ and CoPv02'^'^'^ . 
The physical position of these resistance genes was esti- 
mated to be around 40 Mbp of chromosome 2, based on 
the physical position of the InDel markers. In silico 
analysis (www.phytozome.net) was consistent with the 
location of genes involved in a resistance response in 
thisposition: eleven genes encoding NBS-LRR pro- 
teins, or protein kinases, were annotated in genotype 
G19833 between 39.860 and 40.926 Mbp of chromo- 
some 2 (see Additional file 3 for more detail). 

The location of five resistance genes from Xana on LG 
Pv04 iCo-Sc^'^, Co-3c^-^, Co-Scl'^^^, Co-3c**'''^ and Co- 
all showing a complementary mode of action, 
was also deduced from the genetic analysis. Different an- 
thracnose resistance specificities were mapped in this 
relative position from the bean genotypes Mexico222, 
Widusa, BAT93, JaloEEP558, MDRK, Kaboon, A252 and 
A493 [15,18,19,23,25,26,33]. In genotype Kaboon, one of 
the genes located at this position exhibited a complemen- 
tary inheritance against race 31 [15]. On this position were 



mapped the well-known resistance genes Co-3 and Co-9 
[14]. In this LG, but in a distal position from Co-9, was 
mapped Co-10 from cultivar Ouro Negro as a single gene 
conferring resistance against races 23, 64, 73, 81 and 89 
[22,24]. To date, independence between Co-10 and resist- 
ance genes included in the Co-3 cluster has not been 
clearly established. 

This study confirmed a wide variation in the fungus- 
plant interaction. Up to 22 race-specific genes could be 
involved in the resistance to 11 isolates. Alternative spli- 
cing generates multiple transcript variants from single 
genes and could explain the observed variation in in- 
stances in which recombination among race-specific loci 
was not detected. Alternative splicing was found for the 
RCTl gene in Medicago truncatula that confers resist- 
ance to multiple races of C. trifolli [43]. A complemen- 
tary mode of action indicates the cooperation of two 
genes (or proteins) in the resistant reaction. Combin- 
ation of several proteins in the recognition of the fungus 
can be a mechanism generating variation in the inter- 
action; the number of possible combinations for resistance 
increases. In the guard and decoy models, the effector 
modifies an accessory protein, which may be its virulence 
target (guard model) or a structural mimic of such a target 
(decoy model). The modified accessory protein is recog- 
nized by the NBS-LRR receptor. Under the bait model, 
interaction of an effector with an accessory protein facili- 
tates direct recognition by the NBS-LRR receptor [1,44]. 
Two resistance genes could act in a complementary way 
under any of the three models, with the presence of the 
two proteins necessary to trigger the resistance response. 
Complementary resistance genes required for mediated 
resistance pathways have been described in common bean 
(against the bean rust pathogen Uromyces appendiculatus 
[45]), tomato [46,47], Arabidopsis [48] and barley [49]. 
However, no references were found in a model involving 
three resistance genes which are complementary two-by- 
two. It is important to note that only if the two parents 
differ in the two complementary genes, this type of genetic 
control can be detected. Segregation corresponding to a 
single gene would be expected if the parental lines differ 
in only one of the two complementary genes. In conse- 
quence, the resistance gene detected in a bean genotype 
against a C. lindemuthianum isolate can depend on the 
segregating population investigated. 

Conclusions 

Detailed knowledge of the genetic control of traits (in 
this case, resistance to anthracnose) is relevant for the 
identification of genomic sequences involved in the ex- 
pression using a positional cloning approach, and the 
development of plant breeding programs. Results con- 
firmed that the response in the P. vulgaris- C. linde- 
muthianum interaction was very specific, conditioned by 
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both pathogenic variation of the fungus and the bean 
genotype. Genetic control of the resistant reaction in com- 
mon bean can be controlled by different race-specific re- 
sistance genes depending on the bean genotype. Many of 
these genes were organized as clusters of closely linked 
race-specific genes located in specific chromosome re- 
gions, and these resistance genes can exhibit dominant 
epistatic interaction or complementary mode of action. 
The identification of resistance clusters suggest that most 
of the described alleles for the Co-1 and Co-3 loci are dif- 
ferent and that closely linked loci control the resistance 
reaction against different races or isolates in a bean geno- 
type. All these aspects should be considered in interpret- 
ation of the genetic analysis and the reported results. 

Methods 

Plant material 

A population of 104 Fy RILs developed from the cross 
Xana X Cornell 49242 was used in this study [34,35]. 
The population (XC population) was obtained by single 
seed descent method from individual F2 plants. Xana is a 
bean variety developed at Servicio Regional de Investiga- 
cion y DesarroUo Agroalimentario (SERIDA, Villaviciosa, 
Spain) originated from a cross between the landraces 
Andecha and V203. It is a very large white-seeded line, 
with determinate growth habit and belongs to the fabada 
market class. Cornell 49242 is a very small black-seeded 
line included in the black turtle market class, with inde- 
terminate prostrate growth habit. Cornell 49242 is one 
of the 12 common bean anthracnose differential culti- 
vars used to identify C. lindemuthianum races [12]. The 
remaining 11 differential cultivars were used as a control 
to confirm the identity of the C. lindemuthianum races. 

Inoculation procedure and disease scoring 

Eleven isolates of C. lindemuthianum classified in differ- 
ent races according to Pastor-Corrales [12] were used in 
this study: races 7, 39, 65, 73, 357, 449 and 453, from 
the collection of the Crop and Soil Sciences Department 
(Michigan State University, USA); and races 3, 6, 19 and 
38 from the SERIDA collection. All isolates were ob- 
tained from monosporic cultures maintained in fungus- 
colonized filter paper at -20°C for long-term storage. To 
obtain abundant sporulation, the isolates were grown 
at 21°C in darkness for 10 d in potato dextrose agar 
(DIFCO, Becton Dickinson and Company, Sparks, MD, 
USA). Spore suspensions were prepared by flooding the 
plates with 5 ml of 0.01% Tween 20 (Sigma-Aldrich, St. 
Louis, MO, USA) in sterile distilled water and scraping 
the surface of the culture with a spatula. Inoculations 
were carried out by spraying 8-10-d-old seedlings with 
a spore suspension containing 1.2 x 10^ spores/ml. The 
seedlings were maintained in a climate chamber at 
20-22°C, 95-100% humidity and 12-h photoperiod. 



Responses of the plants were evaluated after 7-9 d using a 
1-9 scale [50]. Seedlings with no visible symptoms (sever- 
ity value 1) or showing limited necrotic lesions (severity 
values 2-3) were considered resistant (R" = resistant reac- 
tion against race X). Seedlings with large sporulating le- 
sions (severity values 4-8) or dead (severity value 9) were 
considered susceptible (S" = susceptible reaction against 
race X). 

The response to a specific race in the XC population 
was evaluated by inoculating all recombinant lines in the 
same test, including 8-10 seedlings per line. The parents 
Xana and Cornell 49242 and the remaining 11 common 
bean anthracnose differential cultivars were also inclu- 
ded in each test. 

Resistance genes were named by using their location 
in anthracnose resistance clusters (Co-cluster), name of 
the isolate or race (in superscript) followed by the bean 
genotype in which the resistance gene was identified. 
For example, a gene conferring resistance to race N in 
Cornell 49242 located in the Co-2 cluster was identified 
as If a resistance gene was not located in re- 

gions in which Co genes were previously mapped, it was 
designated using the LG in which its position was esti- 
mated. For example, a gene conferring resistance to race 
M in Cornell 49242, located in LG Pv09, was tentatively 
named CoPvO^^^. Genes with a complementary mode 
of action were indicated using the letter 'c' after the 
name of the gene (e.g. Co-2c^^'"). 

Genetic linkage map 

A genetic linkage map developed in the XC RIL popula- 
tion was used as a support for the localization of the 
chromosome regions involved in the resistance response 
to different anthracnose races. The map consisted of 
379 loci distributed across 11 LGs [34,35]. MAPMAKER 
Macintosh version 2.0 software [51] was used for the map 
construction using a log of the likelihood ratio (LOD) 
threshold of 3.0 and a recombination fraction of 0.25. 
Marker order was estimated based on multipoint com- 
pare, order and ripple analyses. Distances between ordered 
loci (in centimorgans) were calculated using the Kosambi 
mapping function. The resulting map had 11 LGs, which 
were aligned according to the common bean core linkage 
map using common molecular markers as anchor points. 
LGs were named according to Pedrosa-Harand et al. [52]. 

The XC genetic map included markers tagging the re- 
gions in which anthracnose resistance loci were located: 
markers CV542014 and OFIO530 tagging the anthrac- 
nose resistance cluster Co-1 (LG PvOl) [41]; gene /,/ and 
SW13 marker for the Co-u resistance gene on LG Pv02 
[18,53]; 254-G15F and SW12 tagging the resistance clus- 
ter Co-3 (Pv04); seed protein Phaseolin (Pha) and mar- 
ker SZ4b for the resistance cluster Co-5 (Pv07); marker 
SBB14 for cluster Co-4 (Pv08); and markers SQ4 and 
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SCAreoli which tag anthracnose resistance cluster Co-2 
(Pvll). Recently, a new anthracnose resistance gene 
in common bean, identified as Co-13, was located on 
LG Pv03, linked to RAPD marker OPV2O700 [20]. This 
marker was not polymorphic in the XC RIL population, 
although its relative position on LG Pv03 should be rep- 
resented across the 35 loci forming this LG in the XC 
linkage map. 

Genetic analyses 

Goodness-of-fit of observed to expected ratios was tes- 
ted by chi-square. Additional file 1 summarizes segrega- 
tion ratios expected in a RIL population under different 
hypotheses, considering one to three genes and different 
modes of action. 

To identify the gene(s) involved in the resistance to a 
specific isolate, a systematic genetic analysis was per- 
formed as follows: 

i) Contingency chi-square analyses. Contingency chi- 
square tests of the joint segregation for each scored re- 
sistance with each marker included on the XC linkage 
map were performed. A significant deviation from ran- 
dom segregation would suggest that the chromosome re- 
gion tagged with the marker could be involved in the 
resistance response. Significance thresholds were deter- 
mined using Bonferroni correction from a-level of 0.05 
[54]. First, this analysis was focused on markers tagging 
the main six chromosome regions which included Co 
genes {Co-1 on LG PvOl; Co-2 on LG Pvll; Co-3 on LG 
Pv04, Co-4 on LG Pv08, Co-5 on LG Pv07 and the re- 
gion of LG Pv02 in which the Co-u gene was mapped). 
If a resistance gene was not localized in one of the main 
anthracnose resistance clusters, the contingency chi- 
square test was conducted using the remaining 368 loci 
included in the XC genetic map. 

ii) Direct mapping. When the segregation ratio and 
the contingency chi-square analysis suggested the pres- 
ence of one resistance gene, it was directly included in 
the genetic map. 

iii) Subpopulation analyses. When the segregation ra- 
tio and the contingency chi-square analyses suggested 
the presence of more than one resistance gene, subpo- 
pulation analyses were performed. Two subpopulations 
were established per region, considering the parental ge- 
notypes for the markers that tagged it. Two markers per 
region were used to reduce the possibility of recombin- 
ation events. If a resistance gene was located in the re- 
gion marked, changes in the segregation ratio compared 
to that of the total XC RIL population were expected in 
the established subpopulations. If the region considered 
in the subpopulations was not involved in genetic con- 
trol of the resistance, no change was expected in the seg- 
regation ratio of the resistance, compared to that for the 
complete XC RIL population. 



iv) Genetic dissection. In particular cases, relative posi- 
tion and mode of action of specific resistance genes were 
verified through crosses between selected RILs. These 
lines were selected considering the resistance response 
and their genotypes (Xana or Cornell 49242) for under- 
lying markers located in the putative position of the gen- 
etic map in which the candidate gene was located. 

Additional files 



Additional file 1: Segregation ratios expected in a RIL population 
under different hypothesis. R, resistant; S, susceptible. 

Additional file 2: Observed segregations for resistance to races 65, 
3, 7, 19, 449, and 453 in six subpopulations formed from the XC RIL 
population for three chromosome regions, Co-1, Co-2 and Co-3. 

Additional file 3: Candidate resistance genes annotated from 31 to 
41 IVlbp of chromosome 2 and from 26 to 30 Mbp of chromosome 
9, respectively, using the G19833 genotype sequence available at 
www.phytozome.net. 
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